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We demonstrate that transport spectroscopy of single molecular magnets shows signatures of 
quantum tunneling at low temperatures. We find current and noise oscillations as function of bias 
voltage due to a weak violation of spin selection rules by quantum tunneling processes. The interplay 
with Boltzmann suppression factors leads to fake resonances with temperature-dependent position 
which do not correspond to any charge excitation energy. Furthermore, we find that quantum 
tunneUng can completely suppress transport if the transverse anisotropy has a high symmetry. 
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Introduction. Single molecular magnets (SMM) have 
become famous in the last decade for showing the quan- 
tum tunneling of a single magnetic moment (QTM) on 
a macroscopic scale 1, 2, 3, 4, 5, 6]. These molecules 
are characterized by a large spin 6* > |, a large mag- 
netic anisotropy barrier (MAB) and anisotropy terms 
which allow this spin to tunnel through the barrier. The 
anisotropy is due to spin-orbit effects on the metal-ions 
whose spins couple to form the large magnetic moment. 
Magnetic hysteresis, associated with QTM, was observed 
at temperatures below the MAB 1, 2] for ensembles of 
molecules in a single crystal. Recently, Cornia et al. 
were able to immobilize single SMMs on a gold surface 
through modification of their ligands while preserving the 
magnetic properties of the core. Using this technique 
Heersche et al. Q were able to establish a 3-terminal 
electrical contact and measure the transport through the 
well-known SMM Mnl2, see also 9]. 
In this Letter we show that transport spectroscopy of sin- 
gle molecular magnets can reveal specific features being 
fingerprints of spin quantum tunneling. Even when the 
anisotropy terms which cause QTM have a small effect 
on the energy spectrum they lead to significant changes 
in the non- equilibrium occupations of the magnetic states 
since they allow for a violation of spin-selection rules for 
electron-tunneling. As a consequence, QTM leads to an 
oscillatory behavior of the current and shot-noise with 
increasing bias voltage. Specifically, the interplay of sev- 
eral small rates (quantum tunneling induced rates and 
rates suppressed by Boltzmann factors) leads to nega- 
tive differential conductance and, most strikingly, to the 
occurrence of so-called fake resonances which do not cor- 
respond to any charge excitation energy. The fake res- 
onance's position depends on temperature, and allows a 
clear experimental identification of quantum tunneling 
processes. Furthermore, we show that high symmetry 
(due to the molecular structure) QTM can give rise to a 
complete current suppression. 

Theory. We analyze a minimal model that combines 
the well-known effective spin Hamiltonian description of 
SMMs 'iH m E H IB] with the standard tunneling 
Hamiltonian for the coupling to metallic electrodes. Due 



to the high charging energy it is sufficient to consider only 
two charge states {N = 0, 1) with a magnetic excitation 
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(We employ units h = e = ks = I and energy units 
of meV). For each charge state N the spin has a defi- 
nite value S^^"^ and spin projection \M\ < S^^"^ which 
is maximal in the ground state. The anisotropy terms 
arise due to spin-orbit interaction on the molecule and 
break rotational invariance in spin-space. The lowest 
order easy-axis anisotropy in Eq. Q defines the pre- 
ferred axis in space along which we quantize the spin 
(z-axis). The eigenstates \N,S,M) of Eq. Q have an 
inverted parabolic energy dependence depicted in Fig. 
Cfa). Higher order corrections to the magnetic anisotropy 
barrier are not essential here. It is known experimentally 
15,"li,[ilDj,ll9ll2ai21| and theoretically 14] that the 
anisotropy constants D^^^^ depend on the charge state. 
The transverse anisotropy, Eq. accounts for devia- 
tions from purely axial symmetry. We consider either 
a second or fourth (n = 1 or 2) order term which al- 
lows for tunneling of the spin between states with M 
values differing by 2n. (It is convenient here to deviate 
from the conventional notation E = B2 and C = B4). 
Since a charge-dependent QTM induces only small cor- 
rections in the spectrum of the molecule B2n is taken as 
charge-independent. Transport through SMMs provides 
information on the magnetic structure in more than one 
charge state of the molecule. Therefore we investigate 
the basic possible combinations of the magnetic param- 
eter values for charged SMMs which are scarcely known. 
Also, in a single- molecule junction they may change due 
to mechanical and electrostatic effects. Below we select 
our values from the typical range - 0.01-O.lmeV, 
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10"'^meV and B, - 10""^ - 10"'^meV for 



which magnetic excitations can be resolved in the trans- 
port at electron temperatures below 1 K. Since QTM 



2 



(a) 



N-o (b) 





FIG. 1: (a) Magnetic excitation spectra for two charge states 
with spins S^^^ = 2 > S^^^ =3/2 and D^^^ > D^^\ Since 
typically B2n <^ D'^^^ we label the eigenstates by the approx- 
imately good quantum number M. The dot-dashed line is a 
spin-forbidden transition, all others are spin-allowed, (b) En- 
ergy levels and spin-allowed transitions for S^^^ — S > S'^^^ — 
S — ^. (c) Electron-addition excitations in (V^ate, Vbias) sta- 
bility diagram. Arrows indicate how to construct the diagram 
going along the zig-zag path in (b). Thick/thin lines indicate 
visible/hidden transitions. A transition is hidden when the 
initial state is not yet occupied (by other processes) at the 
transition energy. 



weakly affects the energy spectrum we will label eigen- 
states of H^^^ by the approximately good quantum num- 
ber M, i.e. state \N^S^M) has the largest contribution. 
The electrodes r = L^R are described as elec- 
tron reservoirs with electrochemical potential /i ± 
Vbias/2 and a constant density of states p: H^es = 

Y^rkai^kar - Pr)clf^^Crka' The tunncHng term i^mol-res = 

Yrkja '^j^]a^rka + h.c. dcscribcs charge transfer between 
electrode and molecule (symmetric tunneling barriers). 
Here dj^ adds an electron with spin a to a single-particle 
orbital on the molecule. The coupling to a gate elec- 
trode is included in a shift of the molecular energies, 
such that the charge degeneracy point is at zero bias 
(^ias) and gate voltage (Igate)- For weak tunneling, we 
use a standard master equation approach to calculate 
the non-equilibrium occupations of the molecular states, 
the current and the shot noise 22]. The rates in this 
master equation are calculated in golden rule approxi- 
mation. For the transition S2 si {si being two eigen- 
states of H'^^^^ with energy Ei), we obtain the total rate 
^si,s2 = Yr^lus2 + ^Il7s2 ^^^^ tunneling-in rate 
WIl^s, = '^^pEafriEs^-Es.m^J' andthetunneling- 
out rate I^J;-^ = 2^pE.(l - - ^.J)|Ti;,J^ 

Here, friE) = (e^E-f,,)/T ^ ^yi -g p^^^^- 
tion of reservoir r, T denotes the temperature, and 
^s%2 ~ X]j ^j('^i|^jcrl'^2)- These tunnel matrix elements 
incorporate the spin selection rules and their violation 
for finite QTM. Without QTM, the eigenstates are given 
by l^i) = \Ni, Si, Mi)^ and the tunnel matrix elements 
fulfill obviously the spin selection rule l^*! — 5'2| = 1/2 
and |Mi - M2I = 1/2, in addition to \Ni - A^sl = 1- 
For weak QTM, we decompose the states Si into a 
linear combination of \Ni, Si, M[) states, the one with 
largest contribution being = Mi. Inserting this ex- 
pansion into T!g%2 l^a^s a summation of matrix ele- 



(a) 



(b) 



(c) 



FIG. 2: dl /dVhia.s in gray-scale (gray= zero, white/black = 
positive/negative) as function of V^ate and Vbias- Parameters: 
^(0) ^ 2,5'^^^ = 3/2, D^"^^ = Q.l,D^^^ = 0.01 and T = 0.01. 

(a) No QTM: B2 = 0. (b) QTM B2 = 2 x 10"^ (c) Same as 

(b) except for higher T = 0.015. 



ments with terms (iVi, Si,M[\ tjd]jN2, S2,M^). Us- 
ing the Wigner-Eckart theorem, each of these matrix ele- 
ment can be factorized into an M-dependent Clebsch- 
Gordan (CG) coefficient and a common constant cj. 
Each individual CG-coefficient fulfils l^*! — 5'2| = 1/2 
and \M[ — M2I = 1/2. Obviously the overall spin se- 
lection rule I Ml — M2I = 1/2 can be weakly violated, 
i.e. the corresponding rate is smaller by roughly a factor 
{B2n/ D^^^Y compared to the rates fulfilling the over- 
all spin selection rule. The constants Cj are incorpo- 
rated into a factor F = 2T:p\Y-tjCj\^ common to all 
rates and drop out of the problem, except for setting the 
absolute current and noise scale. We note that, in con- 
trast to customary spin-blockade physics |2^] , a complete 
elimination of the spin projection M from the transport 
problem is not possible due to the MAB and QTM. The 
master equation approach correctly accounts for both the 
non-equilibrium induced by the electron tunneling at fi- 
nite bias voltage and the thermal excitation of molecular 
spin- states. The life-time of the latter is also limited by 
other relaxation processes (spin-phonon interaction, nu- 
clear spins, etc.) which are typically I2J] slower than 
electron tunneling processes (time < 1 ns) and are there- 
fore neglected. Furthermore, the spin-phonon interaction 
may be hindered since the phonon spectrum for a single 
molecule is expected to be less dense than in a bulk sys- 
tem. 

Fake resonances and oscillations. To illustrate the effect 
of the QTM on the transport we first explain the conduc- 
tance map for small spins {S^^^ = 2, S'^^^ =3/2) and con- 
trast the cases B2 = and B2 0. For B2 = the differ- 
ential conductance map is plotted in Fig.I2fa) and we dis- 
cuss the resonance lines running upward. Starting at the 
charge degeneracy point {Vgate = Vbias = 0) and increas- 
ing the bias voltage the current initially sets on due to the 
ground state transitions M = ±2 ^ = ±3/2, see Fig. 
^a). Increasing the bias voltage further brings the tran- 
sition M = ±2 ^ = ±1/2 into the transport energy 
window, without any effect on the current: the rate for 
the process vanishes due to spin selection rules, Wf = 
[dot-dashed in Fig. ^a)]. A reservoir spin-1/2 electron 
can not couple two molecular states with \AM\ > 1/2. 
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Therefore this resonance is hidden [2^ . The current 
only increases when the transition M = ±1 ^ = 
±3/2 becomes energeticahy ahowed. At this resonance 
ah states except A/" = 0, M = become occupied equahy. 
At the third resonance the latter state also becomes ac- 
cessible via M = ^ = ±1/2. In the presence of 
QTM, B2 7^ 0, two additional resonances appear, see Fig. 
EJb), one with positive and one with negative dl/dVhias- 
The appearance of negative dJ/dVbias is related to a slow, 
spin- forbidden transition as follows. For B2 <^ D the 
spin-projection M is only approximately a good quan- 
tum number, i.e. each eigenstate is a linear combination 
of states {|A/', S*, M± 2/c)}/e=o,±i,±2,..., with one coeffi- 
cient {k = 0) close to 1. In the A^ = 1 excited state in 
addition to the state = ±1/2, there is thus a small 
admixture ex ^2/^^^^ of state M = ^3/2. The forbid- 
den transition to the A' = ground states composed 
mostly of M = ±2 [Fig. ffl^a)] is now weakly allowed. 
When it becomes energetically allowed the transition oc- 
curs with rate Wf ~ {B2/ D'^^^)'^T and the current is 
suppressed. This is simply because the occupation of the 
A" = 1 excited states reduces the occupations of the states 
which contribute most to the transport current through 
(fast) spin-allowed transitions. In contrast, the positive 
dl/dVhias line which appears is not related to any addi- 
tion energy of the molecule. The current increase occurs 
when the state causing the negative differential conduc- 
tance (NDC) above is depleted at higher bias via a spin- 
allowed transition M = ±1 ^ ~ (dashed in Fig. 

The rate for this process is Wa ~ Tf{AE — V^^as/^) 
where AE denotes the corresponding transition energy. 
This depletion sets in when in- and out-going rates be- 
come equal i.e. Wa ^ Wf. Due to the small factor in 
Wf this occurs already for Vbias/"^ < where Wa ^ 
rexp(— (A^ — y5ias/2)/T). Equating the rates we obtain 
the resonance condition Vuas/'^. - AE oc Tln(D(i)/52) 
which is substantially shifted from the position expected 
naively {Vbiasf^ — AE). The shift is linear in temper- 
ature and logarithmical in the QTM amplitude |23, |2^ • 
The shift with temperature can be larger than the ther- 
mal smearing as illustrated in Fig. El^c). Thus due to 
the asymmetry between electron tunneling rate constants 
intrinsic to an SMM, transport resonances appear even 
when the molecular level is far away from the electro- 
chemical potential. The strong Coulomb charging effect 
and energy quantization on the molecule are essential to 
this effect since they restrict the transport to sequential 
tunneling through two charge states. 
For SMMs with large spin, ^^^^ = 5^^) ± | > 2, the 
above mechanism leads to oscillations in the transport 
quantities, shown in Fig. [2fa)-(c). At low bias the states 
of the "flatter" A^ = 1 parabola become occupied via 
spin-forbidden transitions, see Fig. |2l(a) and the current 
decreases. The depletion of these states by spin-allowed 
transitions increases the current again. Fig. |31(b). Due 
to the peculiar inverted parabolic energy dependence of 
the magnetic excitations, this sequence is repeated when- 
ever a new A" = excitation state can be occupied, see 
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FIG. 3: Transport oscillations induced by QTM. Parameters: 
5'(o) = 10, 5'^^^ = 9^, D^^^ = 0.1, D^^^ = 0.01, B2 = 2 x 10"^ 
and T = 0.015. (a) d//dVbias for small bias: the A = 1 "flat" 
parabola is mapped out by NDC excitations, (b) dl /dVhias 
for large bias: the A = excitations give rise to positive and 
negative differential conductance, (c) din F/dVbias- fake res- 
onance lines correspond to noise suppression (black) lines and 
terminate at the Coulomb diamond edge, (d) Occupations of 
the = 1 states together with InF and / as function 
of Vbias for V^ate = 0. (c) Currcut suppression due to high- 
symmetry QTM: ^4 = 2 X 10"^, = 0. (f) Spin excitation 
spectrum for the charged state A = 1: the value B^S'^ /D'^^^ 
used in (e) lies beyond the level crossing at approximately 
1. Only the 2nd excited state has a non- negligible admix- 
ture of the maximal M' — ±5'^^-* state which is required for 
transport. 



Fig.^b)-(c). With increasing bias the A" = excitations 
are successively occupied whereas the occupations of the 
A" = 1 excitations, and therefore also the current oscil- 
late^ see Fig. Ol^d). Interestingly, all NDC resonances 
in Fig. ini^a) and (b) correspond to addition energies of 
the SMM (as in our previous example). Most other res- 
onances with positive differential conductance are fake 
since they shift with T and B2. The shot-noise also shows 
oscillations as function of 14ias* in Fig- El (c) the Fano- 
factor F = S/ (21) is plotted. The periodic reductions 
of F occur due to concerted reductions of the noise S 
and simultaneous enhancements of the current I. Their 
positions shift with T and B2 and are fake. The noise 
is super-poissonian, F ^ 1^ due to the presence of slow 
and fast tunnel processes that give rise to large current 
fluctuations [2^. Fig. |31(d) clearly shows that equal oc- 
cupation of all A" = 1 excitations associated with small 
spin-forbidden tunnel constants reduces the current and 
simultaneously leads to stronger fluctuations. Depop- 
ulation of these states by a spin-allowed transition en- 
hances the current and reduces the noise. Importantly, 
in Fig. [3Jc) the (white) lines of enhanced noise persist 
in the Coulomb blockade regime, in contrast to the fake 
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(black) lines of noise reduction. The reason for this ef- 
fect is that the noise in the Coulomb blockade regime will 
only increase when more excitations lie in the transport 
window as shown in ^29] . Hence measuring shot-noise al- 
lows for an identification of fake resonance lines without 
changing the temperature (which may lead to unwanted 
changes in the molecular junction). We note that the 
above is valid also for a weaker magnetic distortion than 
2)(i)//)(o) < 1/2 (which is the requirement for the max- 
imum number of oscillations to occur). 
High- symmetry anisotropy. When the dominant trans- 
verse anisotropy has a high symmetry, i.e. B2/ ^ 
and the ratio BJD^^^ differs for AT = and AT = 1 
com plete current blockade may occur. This is shown in 
Fig.|3(^for ^(0) = 10 and S^^^ = 9^. For BJD^^^ - 
a level crossing occurs between ground an d ex - 
cited state in the charge sector = 1 [Fig. |3(f)| : 
the ground states change from a linear combination of 
{| ^^k)z}k=o,i,... to a superposition of {| T (^*^^^ - 
1) ±4:k)z}k=o,i,..- The latter states have very small tun- 
neling overlap with the = ground state which is a 
superposition of {| ^ 5*^^^ ± "^k) z}k^o,i,... for sufficiently 
small B^/D^^^ <C 1/ S'^. Thus the transport suppression 
at low bias signals a high-symmetry QTM perturbation. 
It can also occur for constant D^^^ ^ D^^^ ^ D when 
B4 changes from smaller than D / S'^ in one charge state 
to larger than this value. However, if the symmetry of 
the QTM is also changed by the charging, i.e. the low 
symmetry QTM becomes important {B2 ^ B/^S'^) in one 



charge state, the current blockade can be lifted, since the 
overlap of ground states is restored. This symmetry low- 
ering may be expected when extra or deficit electrons on 
the SMM are strongly localized on a particular metal ion 
contributing to the total spin. 

Conclusion. Transport spectroscopy of magnetic 
molecules is a challenging task since typically many reso- 
nances are hidden by spin-selection rules or do not corre- 
spond to addition energies of the molecule and shift with 
temperature. Measuring shot-noise allows for an iden- 
tification of misleading excitations without changing the 
temperature. If the total spin values are known, a recon- 
struction of the spectrum from the NDC excitations is 
possible. Then the fake resonances allow the determina- 
tion of the quantum tunneling parameter to logarithmic 
accuracy even though it cannot be resolved directly from 
the thermally broadened excitations. Finally, we showed 
that the transport is even sensitive to the symmetry of 
the magnetic anisotropy of the SMMs. The link we es- 
tablished between transport effects and spin-Hamiltonian 
parameters may be extended down to the microscopic de- 
tails of magnetic molecules with further input from ab- 
initio calculations and energy spectroscopy on charged 
states of SMMs. 
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